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ORING  LOSS  DATA  FOR  SQUARE  LAW 
DETECTORS  FOLLOWED  BY  AN  ORING 
DEVICE  AND  AN  ACCUMULATOR 


INTRODUCTION 

Often  it  is  not  feasible  to  display  all  of  the  data  produced  by  many 
signal  processing  systems.  One  such  signal  processing  system  is  depicted 
in  figure  1. 

In  this  system,  N  displays  are  required  to  process  N  beams  of 
acoustic  data.  In  order  to  reduce  the  amount  of  data  displayed  and 
thereby  reduce  the  amount  of  hardware  needed,  it  is  desired  to  explore  a 
process  that  combines  the  input  data  but  at  the  same  time  results  in  the 
least  signal-to-noise  ratio  (SNR)  loss.  One  process  that  provides  a 
reduction  of  data  is  exclusive  ORing,  a  process  wherein  one  picks  that  single 
channel  with  the  most  energy.  For  the  system  shown  in  figure  1,  exclusive 
ORing  can  take  place  at  either  point  A,  B,  C,  or  0. 

In  earlier  studies  (references  1  and  2),  ORing  at  point  D  was  examined. 
In  the  subsequent  pages,  ORing  at  point  B  is  analyzed.  Two  such  cases  are 
examined:  in  one  case  there  are  N  channels  of  noise,  and  in  the  second  case 
there  are  N  -  1  channels  of  noise  and  one  channel  of  signal  and  noise. 


BEAM  1 


BEAM  2 


BEAM  N 


Figure  1.  Signal  Processing  System 
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An  SNR  is  derived  from  the  statistical  properties  of  each  case.  The  SNR  is 
then  used  to  provide  a  quantitative  description  of  the  effects  of  ORing  at 
point  B. 

Some  related  past  work  on  ORing  at  point  B  is  described  in  reference  3. 
This  analysis  is  based  on  that  work. 


MATHEMATICAL  ANALYSIS 

The  system  of  interest  is  shown  in  figure  2.  The  input  to  the  ORing 
device  is  N  channels  of  acoustic  data,  X,,  X£,  ....  X...  The  N  channels 
contain  either  N  channels  of  noise  or  N  -  1  channels  of  noise  and  one 
channel  of  signal  and  noise.  The  random  variables,  X,,  X2,  ....  X^.,  are 
statistically  independent. 


YES 

NO 


Figure  2.  Generic  Block  Diagram 

The  input  signal-to-noise  ratio  (SNRj)  to  the  ORing  device  in  figure  2 
is  given  by 


SNRj 


=  20  log 


=  20  log  dj  , 


(1) 
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t 


where 

m.  -  m 

dj  *  — — -  =  input  deflection  coefficient 

m^  *  mean  of  the  signal  and  noise  at  the  input  to  the  ORing  device 

mQ  =  mean  of  the  noise  at  the  input  to  the  ORing  device 

a  =  standard  deviation  of  the  noise  at  the  input  to  the  ORing  device. 

There  will  be  an  SNR  loss  when  ORing  N  channels  of  data.  Therefore, 
if  the  same  output  SNR  is  desired  for  N  channels  (N  >  1)  as  for  one  channel, 
additional  SNR  is  needed  at  the  input  to  the  ORing  device.  The  additional 
SNR  needed  is  the  ORing  loss: 


SNR 


Loss 


=  20  log 


(2) 


where 

dT  =  input  deflection  coefficient  for  N  channels  (N  >  1) 
dT  =  input  deflection  coefficient  for  one  channel. 

4 

The  ORing  loss  relative  to  the  detector  input  is  given  by 


SNR 


Loss 


(3) 


The  remainder  of  the  analysis  will  be  devoted  to  the  determination  of 
and  dj^  for  some  fixed  output,  which  will  be  defined  in  terms  of  an  out¬ 
put  deflection  coefficient. 

The  output  of  the  ORing  device  in  figure  2  is  mathematically  defined 
in  the  following  way: 


Y  =  Max  (Xr  X2,  ...,  XN) 


(4) 


Since  one  is  interested  in  the  maximum  of  N  channels,  the  cumulative 
distribution  is  used.  The  cumulative  distribution,  F(y),  for  a  random 
variable,  X,  is  defined  by 
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F(y)  =  probability  (X  <_  y) 


J  f(x)  dx 


(5) 


where 

f(x)  =  probability  density  function 
X  =  random  number 
y  =  real  number. 

The  probability  density  function,  f(y),  is  obtained  by  differentiating 
both  sides  of  equation  (5): 


dF  =  _d_ 
dy  dy 


f ( x)dx  =  f(y) 


(6) 


This  is  a  special  case  of  Leibniz's  rule  for  differentiation  of  an  integral. 

Consider  the  case  of  N  channels  of  noise  at  the  input  to  the  ORing 
device.  The  cumulative  distribution  of  Y  for  this  case  is 


FQ(y)  =  probability  (Y  <_y|  all  noise) 

=  probability  (X.,  X ^ >  •••>  X^  <_y|  all  noise). 

(7) 

Since  the  channels  are  statistically  independent,  equation  (7) 
reduces  to 

FQ(y)  -  P(Xj  <y)  P(X2  <y)  P(X3  <  y)  ...  P(XN  <  y)  .  (8) 

But  P(Xj  <_  y)  =  P(X^  <_y)  =  ...  =  P(XN  <_y).  Therefore, 

F0(y)  =  [P(X  <  y)]N  =  Pjj(y)  .  (9) 

The  probability  density  function,  f  (y),  for  the  noise-only  case  is 
obtained  from  equations  (6)  and  (9):  0 
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f0(y) 


N  1  dP„ 

-  M  17 


■  p0(y) 


Next,  consider  the  case  of  N  -  1  channels  of  noise  and  one  channel  of 
signal  and  noise.  The  cumulative  distribution  of  Y  for  this  case  is 

F^y)  =  probability  (X^,  X2,  ....  XN  <_y|  signal) 

=  P(X2  <  y)  P(X2  <  y)  ...  P(XN  <  y)  .  (11) 

Since  there  are  N  -  1  channels  of  noise  and  one  channel  of  signal  and 
noise,  only  N  -  1  terms  are  equal  this  time.  Therefore, 

Fj(y)  =  [P(Xj  <  y)]N_1  P(X.  <  y) 

=  Po_1(y)  Pj(y)  .  (12) 

where  P0(y)  was  defined  previously  in  equation  (9).  The  probability  density 
function,  f,(y),  for  the  signal  and  noise  case  is  obtained  from  equations  (6) 
and  (12):  L 

dFl 

V*>  *  17 

-  (N  -  1)  P*'2(y)  P  (y)  +  PM_1{y) 

o  dy  VJI  o  dy 


-  (N  -  1)  Pq  (y)  pQ(y)  Pj(y)  +  pj’^y)  Pj(y)  .  (13) 


The  output,  Y,  of  the  ORing  device  is  expressed  in  terms  of  an  output 
deflection  coefficient,  d0Y>  governed  by  equation  (14): 

,  .  MY1  '  UY0 
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where 

UY1  =  mean  va^ue  ^e  S19nal  and  noise  after  ORing 

Uyg  =  mean  value  of  the  noise  after  ORing 

ayg  =  standard  deviation  of  the  noise  after  ORing. 

In  order  to  evaluate  equation  (14),  it  is  necessary  to  compute  uyg,  uyj.  and 
°Y0* 

The  mean  value,  uY,  of  a  random  variable,  Y,  is  defined  in  accordance 
with  equation  (15): 


Uy  =  Y 


L 


yf(y)  dy  , 


(15) 


where  f(y)  =  probability  density  function. 

Utilizing  equations  (10),  (13),  and  (15),  we  find  that  the  mean  values, 
uyg  and  uyj,  for  the  signal  absent  and  signal  present  cases  are 


u YO  =  J yf0(y)  dy  =  N Jy  pj-1(y)  pQ(y)  dy  (signal  absent) 

WY1  •  /yfj(y)  dy  =  J y[(N  -  1)  pj'2(y)  pQ(y)  P^y)  +  pj_1(y)  Px(y) 


(16) 


dy  . 

(signal  present)  (17) 

The  standard  deviation  of  the  noise-only  case,  ayg,  is  defined  in 
accordance  with  equation  (18): 


tO 


3  +>Pyo 


-  IV2  2 

‘  lYo  •  UY0 


(18) 


where 


Uyg  *  mean  value  of  the  noise  after  ORing,  given  by  equation  (16) 


?■/ 


y  fQ(y)  dy  *  N 


/> 


(y)  P0(y)  dy  . 
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The  final  device  in  the  system  of  figure  2  is  the  accumulator.  The 
output  of  the  accumulator,  Z,  is  given  by  equation  (19): 


Z 


(19) 


The  accumulator  output,  Z,  can  also  be  expressed  in  terms  of  an  output 
deflection  coefficient  in  the  following  way*. 


d 


OZ 


(20) 


where 


=  mean  value  of  the  signal  and  noise  at  the  accumulator  output 
UZ0  =  mean  va^ue  t^e  noise  at  the  accumulator  output 
°Z0  =  s  tandard  deviation  of  the  noise  at  the  accumulator  output. 

The  deflection  coefficient  at  the  accumulator  output,  dn7,  can  be 
expressed  in  terms  of  the  deflection  coefficient  at  the  0Ringuaevice  output, 
dgy.  Since  the  random  variables  in  this  analysis  are  independent  and 

uncorrelated,  the  following  is  true: 


UZ0  *  MyY0 
WZ1  *  MuY1 

<4  *  MaY0 

aZ0  =^Y0  ’ 


where  M  -  number  of  terms  in  the  accumulator. 
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Inserting  equation  (21)  into  equation  (20)  gives 


d 


0Z 


PZ1  '  UZ0 
aZ0 


MuY1  '  MuY0  _wff(WYl  ~  UY0\ 

‘  V&YO  V  °Y0  1 

-  >/Md0Y  •  (22) 


The  analysis  of  the  system  in  figure  2  is  now  complete.  A  specific 
example  follows. 


EXAMPLE 

It  is  desired  to  quantify  the  ORing  loss  for  N  =  1  through  8  and  M  =  l, 
8,  16,  32,  64,  128,  and  256  in  figure  2.  In  order  to  compute  the  ORing  loss, 
it  is  necessary  to  determine  the  cumulative  distribution  functions  and  the 
probability  density  functions  in  equations  (16)  and  (17). 

When  no  integration  is  performed  (M  =  1),  the  statistics  used  are  a 
basic  chi-square  distribution  for  the  signal  absent  case  and  a  modified  chi- 
square  distribution  for  the  signal  present  case  (reference  4).  There  are 
two  degrees  of  freedom  for  each  case  and,  therefore,  the  chi-square  distri¬ 
butions  reduce  to  simple  exponentials. 

When  integration  is  performed,  the  input  and  output  statistics  are 
taken  to  be  Gaussian.  Utilizing  Gaussian  statistics,  we  will  now  derive 
the  mathematical  expression  for  the  output  deflection  coefficient,  dg^. 

For  Gaussian  input  statistics,  the  cumulative  distribution  functions 
and  the  probability  density  functions  for  equations  (16)  and  (17)  are 
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Pjx)  *  * 


pp) 


Pj(x)  *  $ 


(23) 


where 


P  (x)  =  cumulative  distribution  of  X  for  the  signal  absent  case 

P^(x)  *  cumulative  distribution  of  X  for  the  signal  present  case 

PQ(x)  =  probability  density  function  of  X  for  the  signal  absent  case 

Pj(x)  *  probability  density  function  of  X  for  the  signal  present  case 

mQ  =  mean  of  the  noise  at  the  input  to  the  ORing  device 

nij  *  mean  of  the  signal  and  noise  at  the  input  to  the  ORing  device 

a  *  standard  deviation  of  the  noise  at  the  input  to  the  ORing 
devi ce 


*(x) 


j  4>(t)  dt 


s  1  /2 

*(t)  3  V27e 

Substituting  equation  (23)  into  equation  (16)  yields,  after  some 
manipulations , 


where 


WY0  3  mo  +  oAN  ’ 

An  ■  N  /*x<(.(x)4N“1(x)  dx 


(24) 
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Substituting  equation  (?3)  into  equation  (17)  yields,  after  some 
manipulations, 


UY1 


=  m^  +  oC^(dj) 


(25) 


where 

CN(dj )  ■  (N  -  1  )J"x<j>(x)$N_2(x)  *(x  -  dj)  dx 

+y*x$N_1(x)  4>(x  -  dj)  dx 
m,  -  m 

j  _  1  o 

dT  -  - 

I  <j 

It  remains  to  determine  aYQ,  the  standard  deviation  of  the  noise-only 
case  after  ORing. 

The  standard  deviation  of  the  noise-only  case,  aYQ,  was  defined 
previously  by  equation  (18).  The  variance,  ag  ,  in  equation  (18)  is  defined 
in  accordance  with  equation  (26): 

2  .  72  2 

°yo  '  Yo  '  uYo 

sJy2  f0(y)  dy  -  uYo  .  (26) 

Employment  of  equations  (10),  (23),  and  (24)  reduces  equation  (26)  to 


vn  =  +  2in  a  A.,  +  a2BM  -  (nr  +  aAj2 


YO 


"0  N 


•  a2(BN  -  A2)  , 


(27) 


where 


vfx 


Bw  ■  N  /  x2  (j)(x)  $N-1(x)  dx 


A^  *  constant  for  a  fixed  N,  defined  previously  in  equation  (24) 
o  *  standard  deviation  of  the  noise  at  the  input  to  the  ORing  device. 
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Therefore,  insertion  of  equation  (27)  into  equation  (18)  yields  the 
standard  deviation  of  the  noise-only  case,  oyg: 

aY0  =  a^BN  "  An)  '  ^28' 

The  deflection  coefficient  at  the  accumulator  output,  dgz»  is  now 
obtained  by  substituting  equation  (24),  (25),  and  (28)  into  equation  (22) 


d0Z  =^d0Y  =  ^ 


yYl  "  kY0 


m  +  oC..(dT)  -  (m  +  cAw) 
-»/n  o  N'  I  o  N 


L(bn  - 


The  statistics  of  the  accumulator  output,  Z,  must  be  evaluated.  Since 
the  input  statistics  were  assumed  to  be  Gaussian  distributed,  the  output,  Z, 
will  also  be  assumed  to  be  Gaussian  distributed.  Therefore,  the  probability 
density  functions  of  Z  for  the  signal  absent  and  signal  present  cases  are 
governed  by  equation  (30): 


(signal  absent) 


rx(z) 


2  -  y, 


(signal  present), 


where  $(z)  3  Gaussian  density  function  defined  previously. 

The  detection  probability,  PQ  and  the  false  alarm  probabili ty,  Pp,  of 
the  system  in  figure  2  is  defined  Tn  accordance  with  equation  (31): 


>D  r,(z)  dz  -  . 

,f  .f  r0(z)  dz  -  .  (*£)  . 
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where 


T  =  threshold  shown  in  figure  2 

r  (z),  r.(z)  =  probability  density  functions  defined  previously  by 
0  equation  (30). 

Consider  the  case  of  T  *  u^.  Por  this  case. 


Pp  *  4>(o)  *  0.5 


*(-doz) 


(32) 


Solving  equation  (32)  for  the  output  deflection  coefficient,  dgZ>  gives 

d0Z  =  "*’1(pF)  for  PD  =  °'5’  (33) 

Therefore,  the  output  deflection  coefficient,  dQZ,  defined  previously 
by  equations  (29)  is  equivalent  to  equation  (33)  for  a  detection  probability 
of  0.5: 


d0Z  =  * 


cN(di)  -  an 


-  A2V 
vbn  \) 


•1 


(Pe)  for  Pr 


0.5. 


(34) 


The  output  deflection  coefficient  for  various  false  alarm  probabilities 
is  listed  in  table  1. 


Table  1.  Output  Deflection  Coefficient  for  Pn  =  0.5 


doz 


10‘2 

2.33 

10"3 

3.09 

1 

o 

f-H 

3.75 

<0 

1 

o 

rH 

4.75 

*This  value  is 

taken  from  reference  5. 

Equation  (34)  can  be  rearranged  in  a  more  convenient  form: 
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W  ‘  AN 


-«f'1(PF) 


for  Pr 


0.5. 


(35) 


This  equation  gives  an  expression  for  the  input  deflection  coefficient 
dj,  required  for  a  specified  Pp  and  PD  *  0.5  as  a  function  of  the  number  of 
channels  ORed,  N,  and  the  number  of  samples  in  the  accumulator,  M.  The  left 
hand  side  of  equation  (35)  can  be  put  in  its  proper  mathematical  form  by 
inserting  the  appropriate  expressions  for  CN( d j )  and  A^  defined  previously: 


W  ‘  AN 


-/«[' 


(N  -  1)  «N_2(x)  4>( x )  $(x  -  dj)  +  4>N_1(x)  <p(x  -  dj) 


N*(x)*N_1(x) 


dx  . 


(36) 


The  constants,  An  and  Bn,  given  in  equation  (35)  are  a  function  only 
of  the  number  of  channels  ORed,  N.  Their  mathematical  expressions,  defined 
previously  in  equations (24)  and  (27),  were  programmed  on  a  PDP-12  computer 
using  OS/8  BASIC.  The  integration  was  performed  using  Simpson's  one-third 
rule  with  automatic  halving  (reference  6).  The  numerical  values  obtained 
for  the  constants,  A^  and  B^,  are  listed  in  table  2. 


Table  2.  Constants  AN  and  B^ 


Number  of 

Channels 

ORed,  N 

an 

bn 

1 

0.00000 

1.00000 

2 

0.56418 

1.00000 

3 

0.84627 

1.27566 

4 

1.02936 

1.55132 

5 

1.16295 

1.80000 

6 

1.26719 

2.02172 

7 

1.35216 

2.22028 

8 

1.42357 

2.39951 

9 

1.48499 

2.56258 

10 

1.53872 

2.71207 

The  input  deflection  coefficient  for  N  *  1,  d^,  can  be  derived 
manually  from  equations  (35)  and  (36)  and  table  2.  The  result  is 


d 


-*_1(PF) 

VR 


(37) 
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where  M  =  number  of  samples  in  the  accumulator. 

The  input  deflection  coefficient  for  N  =  1,  d,  ,  for  a  specified  PF 
and  M  is  listed  in  table  3.  *1 


Table  3.  Input  Deflection  Coefficient  (N  =  1)  for  a  Specified  Pp  and  M 


M 

dT  (PF  =  10'2) 
ll  r 

dr  (P,  =  10'3) 
ll  * 

dT  (PF  =  10"4) 

dT  (PF  *  10"6) 

8 

0.82261 

1.09280 

1.32582 

1.68044 

16 

0.58167 

0.77250 

0.93750 

1.18825 

32 

0.41131 

0.54624 

0.66291 

0.84022 

64 

0.29084 

0.38625 

0.46875 

0.59412 

128 

0.20565 

0.27312 

0.33146 

0.42011 

256 

0.14542 

0.19312 

0.23437 

0.29706 

The  input  deflection  coefficient  for  N  >  1,  di 


is  also  obtained  from 


1TmS 

equations  (35)  and  (36)  and  table  2.  Equation  (36)  was  programmed  on  a 
PDP-12  computer  using  OS/8  BASIC.  The  value  of  the  output  deflection 
coefficient  was  set  equal  to  dg,  =  2.33,  3.09,  3.75,  or  4.75.  The  right- 
hand  side  of  equation  (35)  wasucomputed  manually  for  a  specified  number  of 
channels  ORed,  N,  and  M  =  8,  16,  32,  64,  128,  and  256.  Equation  (36)  was 
then  used  to  obtain  the  required  value  of  the  input  deflection  coefficient 
for  each  N  and  M.  The  integration  in  equation  (36)  was  performed  using 
Simpson's  one-third  rule  with  automatic  halving  (reference  6).  Finally, 
equation  (3)  and  table  3  were  used  to  compute  the  ORing  loss. 


SUMMARY 

The  graphical  results  are  presented  in  figures  3  through  10.  The 
following  conclusions  can  be  drawn  from  an  examination  of  these  figures: 

1.  The  greater  the  number  of  channels  ORed,  the  greater  the  ORing  loss, 

2.  The  lower  the  false  alarm  probability,  the  lower  the  ORing  loss. 

3.  The  more  samples  in  the  accumulator  following  the  ORing  device, 
the  greater  the  ORing  loss. 

4.  Exclusive  ORing  at  point  B  in  figure  1  considerably  decreases 

the  amount  of  data  that  must  be  displayed.  However,  there  will  be 
a  large  SNR  loss  if  integration  is  performed  following  ORing. 

5.  The  overall  system  gain  produced  by  the  integrator  in  figure  1  is 
reduced  by  an  amount  equal  to  the  ORing  loss. 

Finally,  the  analytical  results  for  ORing  at  point  C  (after  inte¬ 
gration)  in  figure  1  are  equivalent  to  the  M*1  curves  in  figures  3  through  6 


SNR  LOSS  (dB)  RELATIVE  TO  DETECTOR  INPUT 
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2  3  4  5  6  7  8 

NUMBER  OF  CHANNELS  ORed,  N 


Figure  4.  SNR  Loss  Due  to  ORing  for  Pn  =  0.5  and  Pc 
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